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Iron-arsenide superconductor Nai-^FeAs is highly reactive with the environment. Due to the 
high mobihty of Na ions, this reaction aflects the entire bulk of the crystals and leads an to ef- 
fective stoichiometry change. Here we use this effect to study the doping evolution of normal and 
superconducting properties of the same single crystals. Controlled reaction with air increases the 
superconducting transition temperature, Tc, from the initial value of 12 K to 27 K as probed by 
transport and magnetic measurements. Similar effects are observed in samples reacted with Apiezon 
N-grease, which slows down the reaction rate and results in more homogeneous samples. In both 
cases the temperature dependent resistivity, pa{T), shows a dramatic change with exposure time. 
In freshly prepared samples, pa{T) reveals clear features at the tetragonal-to-orthorhombic {Ts ~ 
60 K) and antiferromagnetic (Tm=45 K) transitions and superconductivity with onset Tc,ons = 16 K 
and offset rc,o// = 12 K. The exposed samples show T— linear variation of pa{T) above Tc, ons=30 K 
{Tc.off~'26 K), suggesting bulk character of the observed doping evolution and implying the exis- 
tence of a quantum critical point at the optimal doping. The resistivity for different doping levels 
is affected below ~200 K suggesting the existence of a characteristic energy scale that terminates 
the T— linear regime, which could be identified with a pseudogap. 

PACS numbers: 74.70.Dd,72.15.-v,74.25. Jb 



I. INTRODUCTION 

In the generic phase diagram of FeAs-based mate- 
rials, domain of superconductivity has maximum Tc 
close to a quantum critical point of the magnetic or- 
der, see Fig. [l] suggestive of a magnetically mediated 
mechanism of superconductivity.'^^ Quantum fluctua- 
tions lead to a characteristic T— linear temperature de- 
pendence of electrical resistivity, Pa{T), and its evolu- 
tion towards a Fermi-liquid behavior with doping (see 
Ref. ^ for a review). Both these features are observed 
in iron pnictides, most clearly in the electron doped 
Ba(Fei_a;Coa;)2As2 (BaCol22 in the following)!^ and 
in isoelectron doped BaFe2(Asi_a;P2;)2.'^ In hole-doped 
(Bai_^K^)Fe2As2 (BaK122, in the following)P resistiv- 
ity also shows a limited range of close to T— linear varia- 
tion above Tc, however, terminated at high-temperatures 
by a tendency to saturation.'^ Similar slope changes 
are found in bot h in-p lane and inter-plane resistivity of 
self-doped LiFeApESl, and is a dominant feature in the 
intcr-plane resistivity, Pc(T), of transition metal doped 
BaTM12Si4Ui where, by comparison with NMR data^^, 
it was associated with the formation of a pseudogap. Sim- 
ilar association was suggested for the explanation of both 
in- and inter-plane resistivity in BaK122."'^^ 

Nai_5FeAs (Nalll in the following) is isostructural 
to LiFeAs (Lilll in the following); however, in its stoi- 
chiometric form, 5 = 0, it is representative of an under- 
doped part of the superconducting phase diagram, re- 
vealing coexisting superconducting (Tc=12 K), magnetic 



(T™=45 K), and orthorhombic (Ts=55 K) phases .liaHlIl 
In line with this conclusion, superconductivity in Nalll 
can be enhanced with pressure^S and Co-doping.fSS! The 
phenomenology of split structural and magnetic transi- 
tions an d coex isting superconductivity is reminiscent of 
BaCol22.1iSl2ll In contrast with Lilll, Nalll can be pre- 
pared with Na defici ency, allowing for some control of 
the Na stoichiometry.'^SlsSl Due to the high mobility of 
Na atoms at room temperature, Na can be removed from 
the sample by controlled oxidative deintercalation reac- 
tion with the environment .f^Sl Here we use this effect to 
study the evolution of the magnetic susceptibility and 
the temperature-dependent in-plane resistivity, pa{T), in 
the same single crystals of NaFeAs as a function of Na 
content. Our main finding is the observation of T— linear 
resistivity at the optimal doping achieved by the envi- 
ronmental deintercalation, suggesting the universality of 
a quantum critical scenario in iron-pnictide superconduc- 
tors. We find that the temperature range of this T— linear 
dependence is bound from the high-temperature side by 
the pseudogap feature, quite similar to BaK122. The 
pseudogap affects Pa{T) even in the parent composition. 



II. EXPERIMENTAL 

Single crystals of NaFeAs were grown by sealing a mix- 
ture of Na, Fe and As together in Ta tubes and cooki ng 
at 950 C , followed by 5 C/hour cool-down to 900 CW^ 
Samples were stored and transported in sealed containers 
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FIG. 1. (Color online) Generic temperature-composition 
phase diagram of electron-doped iron arsenide superconduc- 
tors, exemplified by Ba(Fei_i:Co2:)2As2.'^ Lines of structural, 
Ts, and magnetic, T^, transitions are split with Co-doping 
and superconductivity has maximum Tc at Xopt ~0.07, close 
to the composition where Ts{x) and T^x extrapolate to zero. 
NaFeAs in its parent state is representative of an underdoped 
part of the phase diagram and the temperatures of its split 
structural and magnetic transitions (circles) correspond to 
Co-doped Bal22 with a;=0.048. With this s-axis shift, ac- 
tual Co-doping phase diagram for Na Fei_iCo2,As (down- 
triangles, data from Ref. [20j ) matches closely the phase di- 
agram of Ba(Fei_i:Co2,)2As2. Up-triangles show Tc and ten- 
tative X position of the samples in which doping is achieved 
by Na deintercalation (this study), working as an electron 
dopant, moves the sample towards optimal doping conditions. 



filled with an inert gas. Preparation of samples for dipper 
tunnel diode resonator (TDR) and resistivity measure- 
ments was done in air as quickly as possible (typically 
about 5 minutes). The leftovers of unreacted sodium 
cause intense reaction on the surface of some samples. 
However, cleaved internal parts of single crystals do not 
show any visible reaction and turns out to be relatively 
stable. After sample preparation, samples were measured 
and stored in a dry, inert environment. 

Samples for in-plane resistivity measurements had typ- 
ical dimensions of (1- 2)x0.5x (0.02-0.1) mm^. All sam- 
ple dimensions were measured with an accuracy of about 
10%. Contacts for four-probe resistivity measurements 
were made by soldering 50 fim silver wires with ultra- 
pure Sn solder, as described in Ref. [5S]. This technique 
produced contact resistance typically in the 10 /if2 range. 
Resistivity measurements were performed in Quantum- 
Design PPMS. 

After initial resistivity measurements, samples were 
left mounted on a PPMS puck in air and measurements 
were repeated after several controlled exposures to air. 
Some samples, after initial preparation, were covered 
with Apiezon N-grease, a technique which was used to 
preserve Lilll in our previous studies.'iSliU In the case of 
Nalll we found that Apiezon N-grease does not prevent 
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FIG. 2. (Color online) DC magnetization of fresh cleaved 
sample measured upon warming after cooling in zero magnetic 
field and applying a 10 Oe field (ZFC-W). Inset: magnetic 
hysteresis loop measured at 5 K in the same sample. 



the sample from losing Na, but slows it down giving a 
better control of sample property variation, especially in 
combination with the ultrasonic treatment. 

DC magnetic measurements were performed with a 
Quantum Design MPMS on a freshly cleaved sample. AC 
magnetic characterizat ion w as performed with a tunnel- 
diode resonator (TDRJE^E^, - a self-oscillating tank cir- 
cuit powered by a properly biased tunnel diode. Sam- 
ples were mounted on a sapphire rod and inserted in 
the inductor coil. The measured frequency shift is pro- 
portional to differential magnetic susceptibility of the 
sample.'^ For quick mounting and measurement proto- 
cols, we used a simplified version of the TDR suscep- 
tometer (a "dipper"), which is inserted directly into the 
transport ^He dewar. The trade-off is reduced stability 
and higher temperature-dependent background as com- 
pared to our high-stability '^He and dilution refrigerator 
versions of the TDR susceptometer. Nevertheless, the 
"dipper" is perfectly suitable to study magnetic signa- 
ture of the superconducting transition, especially when a 
quick turn-around is required. 



III. RESULTS 

Figure [2] shows the temperature - dependent DC mag- 
netization of a freshly cleaved sample measured upon 
warming after cooling in zero magnetic field and apply- 
ing a 10 Oe field (ZFC-W). The magnetic hysteresis loop 
measured at 5 K in the same sample is shown in the inset. 
Both results are quite characteristic of bulk superconduc- 
tivity showing robust screening and significant pinning. 
Moreover, as we show below by direct polarized-light 
imaging and resistivity measurements, this bulk super- 
conductivity coexists with structural and magnetic long- 
range order. 
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FIG. 3. (Color online) Temperature dependence of the nor- 
malized frequency shift in the TDR experiment ( /o is the 
frequency shift at Tc-), for sample of NaFeAs exposed to air 
for the indicated number of hours. Oxidative deintercalation 
of Na increases the onset temperature of the superconducting 
transition from 13 K, for fresh samples, to 24 K for sam- 
ples exposed for 28 hours in air (curves with symbols), only 
slightly broadening the transition. On further exposure, the 
onset Tc rises to almost 28 K, however, additional features 
and significant broadening in the temperature dependence of 
the frequency shift show that the sample becomes inhomoge- 
neous. Inset shows the temperature at half of the transition 
where A///o = 0.5. The shaded area represents proliferation 
of degradation in the exposed samples. 



Figure [3] shows the temperature dependence of the nor- 
malized frequency shift (proportional to magnetic suscep- 
tibility) in single crystal sample of Nalll as a function of 
exposure time to air. Here /o is the frequency shift at T^. 
After initial preparation, the sample shows a quite sharp 
superconducting transition with the onset at Tc '-^IS K. 
Air exposure for up to 28 hours shifts the transition to 
Tc ~24 K, almost parallel to the initial curve and does not 
broaden the transition. This observation suggests that 
the effect is truely bulk in nature, since close to Tc, where 
the London penetration depth diverges, the penetration 
of rf field into the sample is determined by the smaller 
of two characteristic length scales, - skin depth in the 
normal state and characteristic sample dimension, both 
of the order of 100 /im in this case. Upon further expo- 
sure to air, the onset temperature of the superconducting 
transition continues to rise until reaching the maximum 
value of Tc ~28 K, however, the transition broadens, and 
additional features start to appear at lower temperatures, 
suggesting that the sample becomes inhomogeneous and 
starts to decompose into NaFe2As2.'2^ The experiment 
with air exposure was repeated on a total of five samples 
and consistently produced the same result. 

Temperature-dependent resistivity of three "fresh" 
crystals of Nalll is shown in Fig. |4] using a normal- 
ized resistivity scale, p/p{?>QQK). The resistivity value 
at room temperature, /5a(300if), shows notable scatter 
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FIG. 4. (Color online) Temperature dependence of the nor- 
malized in-plane resistivity, pa/pa(300-K') for three samples 
of NaFeAs in a fresh state after initial sample handling 
and contact making. Arrows show features in the temper- 
ature dependence due to magnetic and structural transitions, 
and additional broad slope change features at temperatures 
Ts ^100 K, T2 -180 K and Ti Ri290 K. The inset zooms onto 
the superconducting transition and shows the definitions of 
onset and offset criteria for the superconducting transition 
temperature. 



between 200 to 400 /ificm, presumably due to the ef- 
fect of hidden cracks on the effective sample geomet- 
ric factor .I^^I^ll These values of pa{iQQK) are comparable 
within th e erro r bars to the value of ^300 /ificm, t ypical 
for parenpSEH and hole-doped BaK122 compound^^^Ml^ 
as well as for LilllEil^ and is notably lower than the previ- 
ously reported value of above 4500 pVtcTtiW^ The temper- 
ature dependence of the resistivity of the three fresh sam- 
ples A, B and C is very close, with some difference poten- 
tially coming from the difference in time of air exposure 
during sample preparation and contact making. The air 
exposure can also be responsible for smoother features in 
the temperature-dependent resistivity in this study com- 
pared to the previous study.liSl On cooling, pa (T) of fresh 
samples reveals three broad resistivity slope changes at 
Ti - 290 K, T2 -180 K and T3 - 90 K, a smaU pa{T) up- 
turn on cooling below Ta ~56 K, a sharp downturn below 
Tm ~45 K and a superconducting transition with onset 
at Tc.ons '-^ 14 K (sample A) and offset at Tc.off ~12 K. 
Both onset and offset temperatures were determined as 
crossing points of linear Pa(T) extrapolations above, at, 
and below the transition as shown for Tenons in inset in 
Fig. |4j The onset temperature of resistive transition in 
sample C shifts to 26 K, showing that it is indeed more 
exposed to initial degradation during sample preparation. 

Fig. [5] shows the evolution of the temperature- 
dependent resistivity, Pa(T), in samples of Nalll during 
exposure to air with total exposure time up to 36 hours, 
during which sample properties change homogeneously in 
the TDR experiments, see Fig. [3j The top panel shows 
the evolution of sample resistivity, while the bottom 
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FIG. 5. (Color online) (Top panel) The evolution of 
the temperature-dependent resistivity of the sample of 
Nai_aFeAs during oxidative deintercalation of Na upon ex- 
posure to air. Resistivity of the sample at room tempera- 
ture, pa(3007i'), monotonically increases due to macroscopic 
sample degradation and formation of cracks. Bottom panel 
shows the same data, plotted using a normalized resistivity 
scale, p/p(300isr). The normalization procedure removes vari- 
ation of the effective geometric factor and essentially reveals 
doping - independent resistivity close to room temperature 
and a strong variation below 200 K. The inset in the bottom 
panel zooms onto the superconducting transition range. Ar- 
rows show positions of the special features in pa{T) in the 
parent and 12 hour exposed samples, showing suppression of 
the structural/magnetic transition temperature with doping 
and evolution of the crossover features. 



panel shows the same data plotted using a normalized 
resistivity scale, p / p{?iQQK) , and a zoom of the supercon- 
ducting transition area. It is clear that the environmen- 
tal reaction changes the bulk properties of the sample, 
grossly affecting its normal state Pa{T). On the other 
hand, air exposure monotonically increases pa{iQQK) 
, suggesting development of macroscopic defects (most 
likely cracks) during Na deintercalation. Plotting the 
data on the normalized scale removes the variation of 
the geometric factor due to cracks and reveals very pe- 



culiar transformation of the temperature-dependent re- 
sistivity. First, all four curves coincide above approxi- 
mately 200 K, consistent with the idea that the variation 
of pa{iQQK) value comes predominantly from the varia- 
tion of the geometric factor and that the variation of dop- 
ing level plays only a minor role. Second, below 200 K the 
Pa{T) curve systematically decreases with air exposure. 
Indeed, even unsealed resistivity values decrease caus- 
ing curves crossing in Fig. [5]d. With the exposure time 
of 12 hours, the small resistive increase at Ts becomes 
indistinguishable, while the rapid downturn below 
shifts to lower temperatures as expected for the increased 
doping level, and Tenons rises to 30 K while T^^off rises 
to 20 K. For a 24 hour exposure time, the Pa{T) curve 
shows the sharpest superconducting transition with the 
onset at 30 K, crosses zero at 26 K, and a tail at the 
end of the transition with actual zero reached at 20 K. 
With further exposure increase to 36 hours, the tail be- 
lowthe transition develops further, in broad accordance 
of the characteristic time scale with sample degradation 
in TDR measurements. Fig. [3] 



Together with suppression of the structural transition 
and an increase of the superconducting transition, the 
crossover features in resistivity at T3 become completely 
indistinguishable in a 12 hour curve, while the resistivity 
variation through T2 changes from a slight slope decrease 
on cooling in fresh samples to the onset of a rapid resis- 
tivity decrease in exposed samples. Interestingly, Pa{T) 
becomes T-linear above Tc in samples with a 24 hour ex- 
posure, while the high-temperature end of the T-linear 
range is close to T2. 



In Fig. |6] we show the evolution of the TDR magnetic 
susceptibility in samples stored in Apiezon N-grease and 
subjected to ultrasonic treatment to facilitate Na diffu- 
sion between dipper runs. While grease protects sam- 
ples from air, it does not stop Na deintercalation, which 
now leads to a much higher quality of the superconduct- 
ing transition without additional features and only minor 
broadening. 



Figure [7] shows electrical resistivity of samples covered 
with Apiezon N-grease after initial preparation. Between 
two measurements sample was left on a PPMS puck in air 
for 24 hours. Temperature-dependent sample resistivity 
is shown in the top panel of Fig. [Tj while the same data 
on a normalized resistivity scale, p / p{?iQQK) , and a zoom 
of the superconducting transition are shown in the bot- 
tom panel. Exposure of the sample to Apiezon N-grease 
moderately increases pa{iQQK) and produces a sample of 
much better quality judging by a superconducting transi- 
tion, with Tenons =30 K and Tc,off =26 with actual zero 
resistance achieved at 24 K. 
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FIG. 6. (Color online) Temperature dependence of the nor- 
malized frequency shift in the TDR experiment for a sample 
of Nai-^FeAs stored in Apiezon N-grease and treated in an 
ultrasonic cleaner between successive runs. Oxidative deinter- 
calation of Na increases the onset temperature of the super- 
conducting transition from 13 K for fresh samples to 26 K with 
a slight transition broadening, however, without the appear- 
ance of additional features observed in air-exposed samples, 
Fig.i 

IV. DISCUSSION 
A. Structural/magnetic ordering and resistivity 

Figure [8] shows polarized light images of the freshly 
cleaved surface of the sample used in DC magnetization 
measurements shown in Fig. [2j The top panel was ob- 
tained at 60 K (above the magnetic/structural transi- 
tion), whereas the bottom panel shows an image taken 
at 5 K (well below the transitions) . It reveals a clear pat- 
tern of structural domains, providing direct evidence for 
the occurrence of the tetragonal to orthorhombic struc- 
tural transition in our samples. Considering the up- 
turn anomaly in p{T) at 56 K (structural transition) and 
downturn anomaly at 40 K (magnetic transition) mea- 
sured on another piece of the same crystal, Fig. [7j we 
conclude that freshly cleaved samples exhibit bulk coex- 
istence of superconductivity, magnetism and orthorombic 
structural distortion. 

On the other hand, resistivity in the fresh samples 
rapidly decreases below a temperature of the magnetic 
transition, showing that magnetic fluctuations contribute 
significantly to the inelastic scattering above 



B. Evaluation of the doping level change 

In Nai-A'FeAs, similar to BaFe2As2, doping can be ac- 
complished by substitution of Fe by CcpS!^ P-substitution 
of APl, as well as pressure. Because loss of Na is 
equivalent to Co substitution, both adding extra elec- 
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FIG. 7. (Color online) (Top panel) The evolution of pa{T) of 
a Nai_iFeAs crystal during oxidative deintercalation of Na on 
exposure to Apiezon N-grease. Resistivity of the sample at 
room temperature, pa{300K), increases slightly, presumably 
due to microscopic sample degradation and the formation of 
cracks. The bottom panel shows the same data, plotted with 
normalized resistivity scale, p/ p{300K). The inset in the bot- 
tom panel zooms in on the superconducting transition range, 
showing increase of the superconducting transition tempera- 
ture from Tc,ons~^5 K and rc,o// = 12 K in the fresh state to 
Tc,ons~30 K and Tc,off ~26 K in the grease exposed state. 



tron, it is natural to expect that similar compositional 
variations are involved. In Fig. [TJwe compare the phase 
diagram of BaCol22 with that of NaFeAs, using the 
structural transition temperature as a matching crite- 
rion. This way Nai^^FeAs in its "fresh" state corre- 
sponds to a;=0.048 in BaCol22, while at optimal doping 
Nai-^FeAs corresponds to a;=0.07 in BaCol22. This dif- 
ference, Ax=0.022 is very close to the actual Co-doping 
required to induce the highest Tc=25 K in Nai-^FeAs.'^Sl 
This comparison suggest that S in our samples is of the 
same order, S sa0.02. 
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FIG. 8. Polarized light images of the fresh cleaved surface 
of a Nai-^FeAs sample, taken above structural transition at 
60 K (top panel) and at 5 K, the base temperature of our 
experiment (bottom panel). The image shows a pattern of 
structural domains, with domain walls running along the [100] 
tetragonal direction. 



C. Evolution of the temperature-dependent 
resistivity in 111 compounds 

Figure [9] shows temperature-dependent resistivity in a 
111 system. We plot the data for fresh NaFeAs, grease 
treated optimally doped Nai_aFeAs and LiFeAs, repre- 
sentative of the overdoped compositions. For compari- 
son, we show temperature dependence in optimally doped 
BaCol22. Note that despite macroscopic damage of the 
crystal during deintercalation, non-stochiometric compo- 
sition shows quite a high residual resistivity ratio, with 
p{300K)/p{Tc) w6 and p(300X)/p(0) w25. This is no- 
tably higher than in BaCol22 with p{300K) / p{Tc) «3, 
suggestive that doping into the Fe sites introduces much 
stronger disorder than the formation of Na vacancies. 
As can be seen from the inset in Fig. |9) at the opti- 
mal doping, p(T) in Nalll is purely T— linear above T^, 
which should be contrasted with a quad ratic tempera- 
ture dependence observed in LiFeAs.l^Hl] In the case of 
Nalll the temperature range of T-linear dependence is 
bound from the high-temperature side by a downward 
slope change on heating, similar to the case of a hole 
doped BaK122. 

Our findings of doping-independent resistivity at high 
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FIG. 9. (Color online) The evolution of the temperature- 
dependent resistivity, plotted on a normalized resistivity scale 
p/ p{300K), in a 111 system. The data for parent NaFeAs 
(curve 1), grease treated Nai-^FeAs (curve 2) and LiFeAs 
(curve 3) are representative of the underdoped, optimally 
doped and overdoped regimes, respectively. The inset shows 
a zoom of the low-temperature range for curves 2 and 3, 
the line is a linear fit through p{T) of grease-deintercalated 
NaFeAs. For reference we show the in-plane (curve 4) and 
inter-plane (curve 5) temperature-dependent resistivity of op- 
timally doped Ba(Fei_i:Coa;)2 As2, a;=0.074, showing T-linear 
resistivity at low temperatures and a pseudogap crossover for 

Pc(rjSi. 



temperatures, which are particularly clear for air treated 
samples, are very difficult to explain in a model of ad- 
ditive contributions of two different bands into conduc- 
tivity. We suggest that a broad crossover in p{T) at 
around 200 K has a similar origin as a broad crossover 
in c-axis transport in BaCol22, where this effect corre- 
lates well with the NMR observations of a pseudogap and 
with the domain of T— linear magnetic susceptibility.^ 
It is interesting that the pseudogap crossover does not 
affect in-plane transport in electron-doped BaCol22, 
but strongly affects electron-doped Nai-^FeAs and hole 
doped BaK122. Finally, we point out a very pronounced 
branching in the temperature-dependent resistivity, coin- 
ciding with the end of the pseudogap resistivity crossover. 
This feature suggests that the electronic structure of the 
compounds is affected at temperatures much higher than 
the temperatures of structural and magnetic transitions, 
a feature hard to reconcile with simple spin density wave 
ordering models. It brings some similarity to a pro- 
nounced strain-induced anisotropy well above structural 
and magnetic transition^^H jjj ^j^g equivalent doping range 
in BaCol22. 



V. CONCLUSIONS 

We used the naturally occurring oxidative deinter- 
calation procedure to fine - tune the doping level of 
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the Nai-^FeAs system from underdoped to optimally 
doped. Analyzing DC and AC magnetization, direct op- 
tical imaging and temperature - dependent resistivity, we 
conclude that freshly cleaved samples exhibit bulk coexis- 
tence of superconductivity, magnetism and orthorombic 
structural distortion and can be placed on the under- 
doped side of the effective T — x phase diagram. The 
oxidative deintercalation results in a shift of the com- 
position towards the effective optimal doping. Measure- 
ments of the in-plane resistivity show that the difference 
between fresh {5 = 0) and non-stoichiometric states of 
the crystals starts at about 200 K, much higher than 
the temperature of the structural/magnetic transitions, 
Ts ^ 56 K and Tm ^ 40 K. This temperature is of 
the same order as the temperature of a broad crossover, 
found in resistivity measurements in all 111 compounds. 



which we relate to the formation of a pseudogap. At the 
optimal doping, the resistivity of Nai-^FeAs shows an 
extended range of T— linear behavior, in line with the 
expectations of quantum-critical scenariJS, however, this 
range is bound on the high-temperature side by the onset 
of a pseudogap crossover. 
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